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Abstract. The possibility to determine the pA scattering length from the final-state interaction in the 
reaction pp — > pK + A is investigated experimentally. From a transversely polarized measurement, the K + 
analyzing power (^4jv) which, in principle, allows one to extract the spin triplet scattering length is studied. 
An unexpected energy dependence of the forward/backward symmetric part of is found. The influence 
of iV*-resonances on the pA invariant mass spectrum is investigated by exploiting the large acceptance for 
the process pp — > pK + A — > pK + piv~ and is found to be the main source of uncertainty for determining the 
pA scattering length. 

PACS. 13.75.-n Hadron-induced low- and intermediate-energy reactions and scattering (energy < 10 GeV) 
- 13.75.Ev Hyperon-nucleon interactions - 25.40.Ve Other reactions above meson production thresholds 
(energies > 400 MeV) 



1 Introduction 

While nucleon-nucleon scattering has been precisely mea- 
sured and accurately described up to beam kinetic energies 
of 3 GeV [l] , the situation is much worse for hyperon- 
nucleon scattering. Recently, a dispersion relation tech- 
nique was developed by Gasparyan et. al. [2] that allows 
one to extract the pA scattering length from final-state 
interactions in processes with high momentum transfer like 
pp —¥ pK + A. It requires data where the pA system is in a 
specific total spin (S=0 or S=l). In an experiment with 
transversely polarized beam, the symmetric component of 
the K + analyzing power (Am) can be used to disentangle 
the spin triplet component of the pA interaction. Then, 
the scattering length can be extracted from the shape of 
the corresponding pA invariant mass spectrum. 

*-resonances were found to have a significant influence 
on the pKA production cross sections [3j[4] . In Ref. [2] it is 
pointed out, that a Dalitz plot analysis should be performed 
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to check whether the area of final-state interactions is 
overlapping with those resonance structures. 

Due to its unique feature of 47r acceptance for the final 
state, the COSY-TOF spectrometer is ideally suited to 
measure the pp — > pK + A reaction and provide the data 
needed for the extraction of the scattering length and the 
Dalitz plot analysis. It has recently been upgraded for 
improved event reconstruction capabilities. These should 
allow to achieve an experimental precision that competes 
with the accuracy of the extraction method. Additionally, 
the polarized proton beam from the COSY accelerator 
gives access to the polarization observable An. 

In this paper we apply the extraction method to data 
with a beam momentum of 2.95GeV/c, corresponding to 
an excess energy of 203.7 MeV. Systematic effects due to 
the influence of N* -resonances are quantified by separately 
analyzing different regions of the Dalitz plot. We observe 
an unexpected energy dependence of the K + analyzing 
power, and consequently it is not possible to determine the 
spin triplet scattering length at this beam momentum, with 
the amount of data available. Nevertheless, this finding 
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Fig. 1. Schematic view of the COSY-TOF detector, including 
the new Straw Tube Tracker (STT) and Silicon Quirl Telescope 
(SQT). 



is interesting and further experimental and theoretical 
studies are necessary. 



2 Method 

2.1 Experimental Setup 

The new experimental setup of COSY-TOF is sketched 
in Fig. [l] Segmented scintillators close to the target and 
at the inner side of the main vacuum vessel are used for 
triggering, time of flight and dE/dx measurements. New 
detection subsystems have been installed to improve the 
event reconstruction precision and efficiency: the Silicon 
Quirl (SQT) close to the target and the Straw Tube Tracker 
(STT) [5] inside the main vacuum vessel. The STT has 
nearly Air acceptance for the pp — » pK + A — > pK + pir pro- 
cess and detects all charged particles in the final state, 
thus it is the most important subsystem for this analysis. 
It consists of 2704 individual straw tubes combined into 
13 double layers normal to the beam axis with 3 azimuthal 
orientations. The spatial resolution achieved under experi- 
mental conditions of the individual straws has been shown 
to be a « 150 um with an efficiency better than sa 98%. 
For the pK + A final state this results in a m « 1.1 MeV/c 2 
resolution for the pA invariant mass at a reconstruction 
efficiency times acceptance of 25% [6] for the charged decay 
mode. 



2.2 Determination of the Spin Triplet Scattering 
Length (a t ) 

In Ref. [2] a dispersion relation is derived that connects the 
pA scattering length to an integral over the pA invariant 
mass spectrum. The spectrum is integrated over the pA 
invariant mass (m p A) from threshold (mo = m p + ttia) 
up to m max = m + 40 MeV/c 2 , because the pA system is 
required to be in an S-wave. From simulations it is argued 
that, with this limit, the accuracy of the method is 0.3 fm. 



The application of the method requires observables 
where the pA system is in a specific total spin state. In 
Appendix B of Ref. [2] it is pointed out that from a mea- 
surement with a transversely polarized beam observables 
are accessible to which only the spin triplet part of the pro- 
duction amplitude contributes. Therefore, the spin triplet 
scattering length a t can be determined. To show that, the 
K + analyzing power is expanded in terms of the associated 
Legendre Polynomials of degree I and order m: 

^(cos^mpyi) = a(m pA )P^ (cos 9^)+ /3(m pA )P 2 (cos6£). 

(1) 

Here, 6^ is the polar angle of the K + in the center-of-mass 
system and the function P* (P 2 ) is forward-backward 
symmetric (asymmetric). It turns out that higher order 
contributions to Eq. ([I]) are negligible in the current exper- 
iment. This implies that the coefficient a results from the 
interference of amplitudes that correspond to the K + being 
in an s- and p-wave, respectively, while the coefficient f3 is 
due to an interference of the s- and d-wave amplitudes [7j. 

Under the assumption that the pA system is in an 
S-wave, only spin triplet amplitudes contribute to a (see 
Appendix B of [2] for details). This follows from the fact 
that different spin states in the final state do not interfere 
and that, in addition, for a spin zero pA final state, only 
even K + partial waves are possible^] Consequently, there 
is no contribution of the pA singlet state to the K + p-wave, 
hence to a. Therefore, the integral 
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contains only contributions from pA spin triplet states 
because the forward-backward asymmetric /3-term cancels 
out. Here, p pA denotes the relative momentum between 
p and A and d 2 a/dcos9^dm pA is the double differential 
production cross section. In Eq. ^ the full (spin-averaged 
and angle-integrated) production amplitude, 



\M(m pA )f 



1 da 

PpA dm pA 



(3) 



was introduced, which appears since da /d cos 9 is isotropic. 
The latter is consistent with the exclusion of higher partial 
waves in Eq. ([!]). 

Using the parametrization: 



t(m pA )\\M(m pA )\ 2 
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(4) 



For a pA system in an S-wave and spin zero state one has 
for the total angular momentum J tot = Ik, with Ik for the 
angular momentum of the kaon with respect to the Ap system. 
The parity of the final state thus reads n; = (— 1) 1k+1 . Because 
the parity of the initial state is given by the relative angular 
momentum Li of the protons as 717 = (— l) Li , it follows from 
parity conservation that L t 7^ J tot . Therefore, the initial state 
can only be in spin 1 state. Because then the Pauli principle 



demands that (- 
even. 



-1, it follows from Li 7^ 



that 
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with the free parameters C'i, the spin triplet scattering 
length a t can be obtained from: 



o t (C 1 ,C7 2 ) = -^C? 



(™L - m o) 



(5) 



m p m A J (m^ ax - Cf )(m§ - C|) 3 ' 



The independence of the scattering length from Co reflects 
the fact that only the shape of the final-state interaction 
plays a role — this is why a proportionality of the spin 
triplet scattering amplitude to |a(m py i)| \M(m p A)\ 2 is suf- 
ficient. 



2.3 Determination of the K + Analyzing Power A 



N 



The analyzing power Am is a measure for the left/right 
asymmetry 6lr of the K + differential cross section and is 
defined as: 



P 



cos{4>) ■ A N {6* K Y 



(6) 



with the beam polarization P. The asymmetry is deter- 
mined from 



L(pZ,<f>)-R(% 



K ' 



L{6*,4>) + R{0* K , 

7T 7T 
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(7) 



where 



L{e* K ,4>) 

and R{6^,4>) 



^/N+(0)-N-(0 + n) 



(8) 



Here, N ± ((j)) is the number of events with spin up (+) 
and spin down(— ) projectiles at the azimuthal angle <fr. By 
multiplying the number of events on the opposite sides of 
the detector and opposite spin states, systematic effects 
from asymmetries in detector acceptance are canceled to 
first order. The data has been divided into eight bins in 
the K + azimuthal angle <f>. With the known analyzing 
power and the measured asymmetry in pp — > pp elastic 
scattering the beam polarization was determined to be 
(61.0 ± 1.7) %. For that the pp analyzing power was taken 
from the SAID partial wave analysis [T] . The polar angular 
dependence is in good agreement with SAID and with a 
previous measurement by EDDA [8j. 



2.4 Data Analysis 

The primary signature of a good event is the combination 
of two primary tracks from the target and two tracks from 
the delayed weak decay of the A particle. Additionally, the 
A decay plane contains the primary vertex position. After 
selecting events that fulfill these criteria, a kinematic fit is 
performed which minimizes the Xkm with respect to the 
measured track to wire distances in the STT. To reject 
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Fig. 2. Missing energy spectrum at the primary vertex. From 
all events with a successful kinematic fit (dashed line) a sub- 
sample is selected (blue solid line) as described in the text. Also 
the complementary sub-sample is shown (red dotted line). 



background processes a threshold is set on the reduced 
chi-square: 

XL/NDF < 5 (9) 

Additionally, events with a minimum distance sa between 
the production and decay of the A are selected 



sa > 3 cm. 



(10) 



Furthermore, the laboratory angle between the A and its 
decay proton is required to fulfill the condition: 



/.(A,p) > 3mrad, 



(11) 



in order to reduce instrumental background from events 
with multiple primary tracks. 

The effectiveness of these selection criteria is evaluated 
by the distribution of the missing energy (-E^iss) at the 
primary vertex with respect to a pK + A final-state hypothe- 
sis, before the kinematic fit. This is shown with the dashed 
line in Fig. [2] for all fit events with s A > 1 cm. The sig- 
nal peak around MeV lies on top of a continuum from 
the instrumental background. By applying the criteria in 
Eqs. (9pT ) only events in the peak are selected (solid). At 
Pbeam = 2.95 GeV/c a total sample of about 42, 000 events 
is obtained. 

Physical background remains from the process 
pp — > pK + U° — > pK + A^y, where the unmeasured 7 carries 
away w 77 MeV. This results in a deflection of the A by 
pa 2° in the laboratory frame. The event topologies are 
therefore similar. Studies of MC generated events have 
shown that the contamination of the event sample under 
the conditions (9pT) is < 5%. Therefore, it is neglected in 
the following analysis. 



3 Results 
3.1 Dalitz Plot 

The Dalitz plot of the selected event sample is shown 
in Fig. [3] It is corrected for detector acceptance with 
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Fig. 3. The Dalitz plot of the reaction. In addition to a promi- 
nent enhancement in rr? vA at the NE threshold the FSI inter- 
action with influences of iV*-resonances is visible in the upper 
left corner. 



MC generated events. The complete kinematic acceptance 
of the COSY-TOF detector is evident. The Dalitz plot 

at the NE threshold 



10 



and 



density is strongly enhanced in m 2 A 

at 4.53GeV 2 /c 4 . This has been observed before [9 
could be interpreted as a NE-pA coupled channeTeffect. 
The high resolution available in this measurement makes 
an analysis of the shape, position and strength of this 
structure interesting; however that is beyond the scope of 
this report. It is analyzed in more detail in Ref. 

The pA final-state interaction can be seen as an en- 
hancement at small nip a- The increasing differential cross 
section for decreasing ttika (see Fig. [3 ) can be explained 
by the influence of the resonances iV(1710) and/or iV(1720) 
I3jp]. In the Dalitz plot these are located around m 2 KA « 
2.93GeV 2 /c 4 . However, due to their width of more than 
100 MeV they do not appear as narrow structures. 



3.2 Effective Scattering Length 

In Fig. [4] the pA invariant mass spectrum is shown. Since 
the time integrated luminosity of the event sample is not 
needed for this analysis, only the number of measured 
events (N) scaled with the detector acceptance times re- 
construction efficiency (A) is given. The quantity A has 
been determined with Monte Carlo studies and is shown 
on the bottom of the figure. It is noteworthy, that the 
detector acceptance is nearly constant over a wide m p A 
range but varies between 27 % and 10% close to threshold. 

For comparison, an arbitrarily scaled three-body s- 
wave phase space distribution is shown with a solid line. 
The NE threshold positions of m nE + = 2128.9 MeV/c 2 and 
Trips — 2130.9 MeV/ c 2 are marked with the vertical lines, 
respectively. In this region a strong enhancement appears. 
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Fig. 4. The spectrum of m p A corrected for acceptance (A) as 
it is given on the bottom. An enhancement is visible at the NE 
thresholds (vertical lines). An arbitrarily scaled phase space 
distribution (solid line) is shown to guide the eye. The red solid 
line is a fit to the data as described in the text. 



There is also an enhancement of the production cross 
section visible at small m p A due to pA final-state interac- 
tions. To test the method for the extraction of the scat- 
tering length the differential cross section, calculated from 
Eq. Q, convoluted with the detector resolution is fit to 
the data in the range (0, 50) MeV/c 2 above threshold with 
a x 2 /NDF = 0.32 (solid line). 

From a likelihood analysis [2 f 6j of the highly correlated 
parameters we obtain a — (—1.25 ± 0.08 ± 0.3) fm. Here, 
in the second and third term the uncertainties from statis- 
tics and of the theoretical method are given, respectively. 
Because the incoherent sum of spin singlet and spin triplet 
pA scattering enters M(m p A) and the relative weights are 
unknown, a is referred to as effective scattering length. 
Note, that it is not a spin average. 

A variation of the upper limit of the fit range between 
(40, 50) MeV/c 2 yields a stable result within the statistical 
error. At an upper limit of 60 MeV/ c 2 the absolute value 
of the effective scattering length increases by « 0.3 fm. 
Although this is still within the systematic error, it might 
be connected to an increased importance of higher partial 
waves in the {pA} system or to influences of the enhance- 
ment around the NE threshold. 

To quantify the influence of A^*-resonances on the mea- 
surement, we apply the same method as above to two 
separate m^/i-regions of the Dalitz plot. For a better com- 
parability Fig. [5] presents the ratio of the spectrum to 
the arbitrarily scaled phase space (circles), which is pro- 
portional to |Af(m p/ i)J 2 . Therefore, a fit of the exponen- 
tial function of Eq. Q is shown (solid line). The data 
are normalized to have an average \M(rn p A)\ 2 — 1 in 
the range (2090 < m pA < 2110) MeV/c 2 . \M(m pA )\ 2 is also 
given for (2.590 < m 2 KA < 3.176) GeV 2 /c 4 (triangles down) 
and (3.176 < m 2 KA < 3.287) GeV 2 /c 4 (triangles up). These 
two regions are referred to as the lower (£r) and upper (ur) 
ranges, respectively. Their boundaries are chosen such that 
both regions include m pA = m p + m A . To improve read- 
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Fig. 5. The quantity |M| 2 (circles) fit with the exponential 
function of Eq. Q (solid line). It is also shown for the upper 
(triangles up) and lower region (triangles down) of the Dalitz 
plot as described in the text. For better readability the latter 
are shifted +1.5 units on the y-axis, respectively. 



ability the data are shifted by +1.5 units along the y-axis. 
It can be seen that the measured strength and shape of the 
final-state interaction varies significantly between all three 
samples. By applying the same method as above, for the 
lower range we obtain a lr = (—0.86 ± 0.06 ± 0.3) fm and 
a ur = (—2.06 ± 0.16 ± 0.3) fm for the upper range, respec- 
tively. The difference of 1.20 fm shows that the systematic 
effect of iV*-resonances dominates the precision of the de- 
termination of a. This result agrees with our observations 
for the Dalitz plot in Fig. [3j At low m p A the density clearly 
deviates from a homogeneous distribution along rn 2 KA , lead- 
ing to a tilted shape of the enhancement from final-state 
interactions. 



3.3 K + Analyzing Power 

In Fig. [6] the analyzing power of the K + is shown for the 
whole TOpyi range. The parameters a and j3 of Eq. (fll) are 
fit to the data (solid line). The good quality of the fit, 
X 2 /NDF = 0.54, justifies the exclusion of higher order 
contributions and therefore supports the assumption of an 
isotropic da jd cos 8* . The latter was also verified from a di- 
rect measurement of this distribution [4j[6] . The symmetric 
(dotted line) and asymmetric contributions (dash-dotted 
line) are also shown separately. 

In Fig. [7] the event sample is binned in m p A and the fit 
results for — a(m p A) (circles) and j3{rn p A) (triangles) are 
given, respectively. The asymmetric part of An approaches 
zero for higher m pJ \. This is reasonable because in this 
region the K + has the lowest momentum which could be 
insufficient for s- and d-wave interference. The symmetric 
part on the other hand is close to 25% on the high end of the 
spectrum and reduces to < 11 % with 3cr confidence on the 
low end. This means that with the present quantity of data 
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Fig. 6. The K + analyzing power for the whole m p A range. The 
fit (solid line) includes the associated Legendre Polynomials 
(Pi) (dotted line) and (Pf) (dash-dotted line). 
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Fig. 7. The coefficients of the forwards-backward symmetric 
(circles) and asymmetric (triangles) contribution to the K + ana- 
lyzing power as a function of m p A ■ To high m p A the asymmetric 
contribution vanishes because the K + momentum is too low 
for s- and d-wave interference. Surprisingly the symmetric part 
vanishes for low m v A- 



the dependence of a on m pA cannot be determined with 
sufficient precision to extract the spin triplet scattering 
length with Eq. 



4 Discussion 

The observed behavior of the analyzing power is unex- 
pected because in the reaction pd — > d?r + , that has only 
spin 1 in the final state and consequently the same selection 
rules, a value around 25 % for the symmetric contribution 
has been observed 
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Our finding is compatible with the 
hypothetic al ab sence of spin triplet scattering. As discussed 
in Section 2.2 in that case the even K + orbital angular 
momentum contributions would be forbidden at small m p A- 
Then it follows that d would be the 1 Sq scattering length 
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a s . Current models 13 - 15 predict (—1.2 < a t < — 1.8) fin 
and (—1.8 < a s < — 2.6)Tm. Their parameters are fit to 



and tested to predict a suf- 
in the Hypertriton 17 18 



pA — > pA elastic scattering [16 
ficient binding energy of the / 
Since the iV(1710) and Y(1720) lie below the lower range 
in ttika, one can assume a weaker effect in the upper range. 
Therefore, if we interpret the extracted a ur as a lower 
boundary for the spin singlet scattering length, it would 
fit into the theoretical picture. 



This would also agree with a result from HIRES 19 
that finds a = — 2.41^25 at Pbeam = 2.7GeV/c. In the com- 
plete absence of spin triplet scattering the results would be 
comparable. Also, note that in Ref. [20] it is demonstrated 
that the method applied in Ref. 



19 overestimates the scat- 



tering length by approximately 0.4 fm. Moreover, the effect 
of 7V*-resonances on their result is unclear, even though 
earlier measurement of COSY-TOF indicate a reduced 
influence on the Dalitz plot [3] at that Pboam- 

The behavior of the analyzing power can also be ex- 
plained in the presence of spin triplet pA scattering: As 
only imaginary parts enter the analyzing power, a phase 
cancellation is possible. The viability of this hypothesis 
strongly depends on the upper limit one can set for the 
symmetric part of the analyzing power. It is therefore im- 
portant to collect higher statistics and measure at different 
beam energies. It should also be investigated how other 
polarization observables, e.g. the A polarization, can be 
employed to put quantitative constraints on the presence 
of spin triplet pA scattering. 



5 Conclusion 

The effective pA scattering length has been determined 
with an experimental precision that surpasses the uncer- 
tainty of the extraction method. A test of the influence 
of Y*-resonances revealed that they impose the largest 
uncertainty on the result. Therefore, to determine the pA 
scattering length in the reaction pp — > pK + A these influ- 
ences require further study. For that, data from different 
beam momenta is highly desirable. 

The K + analyzing power has been measured as a func- 
tion of the pA invariant mass. The vanishing symmetric 
contribution to the analyzing power at low values of m v A 
prohibits the extraction of the spin triplet scattering length 
with the present quantity of data. The hypothesis of exclu- 
sive spin singlet scattering was discussed as an explanation 
for this unexpected behavior. For a decisive study, higher 
statistics are needed. Especially, the possibilities to exploit 
the A polarization should be investigated. 



We would like to thank the COSY operation crew for 
providing excellent quality beams. 
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